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Metal-insulator transitions in the paramagnetic phase of the two dimensional square lattice Hub- 
bard model are studied using the dynamical cluster approximation with eight momentum cells. 
We show that both the interaction-driven and the doping-driven transition are multi-stage and 
momentum-sector specific, with Fermi liquid metal and fully gapped insulator phases separated by 
an intermediate phase in which some regions of the Brillouin zone are gapped while others sustain 
gapless quasiparticles. We argue that this is the coarse-grained version of a gradually shrinking arc 
or pocket. A pronounced particle-hole asymmetry is found. 

PACS numbers: 71.10.Fd, 74.72.-h, 71.27.+a, 71.30.+h 



The metal-insulator transition in strongly interacting 
Fermion systems is one of the central issues in condensed 
matter physics The transition observed in the high- 
T c cuprate family of materials is particularly interesting. 
As the doping S (here defined as the difference of carrier 
concentration from 1 carrier per Cu02 unit) is reduced 
to zero the materials become less metallic, with an in- 
sulating phase occurring at zero doping. Experiments, 
including analyses of the temperature dependence of the 
London penetration depth 0, H|, normal state Q and 
superconducting Q Raman spectra, c-axis conductiv- 
ities 0], photoemission [§| and scanning tunneling mi- 
croscopy Q (for a review see (l(|) suggest that the tran- 
sition is driven by a doping-dependent decrease in the 
size of the region of momentum space which can support 
well-defined mobile quasiparticles. Understanding this 
behavior has proven theoretically challen ging . Weak cou- 



11, 12] indicate 



pling renormalization group calculations 
a parameter regime in which interactions involving elec- 
trons near the (±7r/2, ±7r/2) points flow to weak coupling 
while interactions important for electrons near the (ir, 0) 
and (0, 7r) portions of the Fermi surface flow to str ong 
coupling for reasons relating to antiferromagnetism [ll| 
or Mott physics [l2j]. However, the basic result is a flow 
to a strong coupling regime where the equations no longer 
apply and the physics must be studied by other methods. 

Dynamical mean field theory (DMFT) [HEil 

has pro- 
vided valuable insights into the physics of strongly cor- 
related materials. Its original single-site version averages 
physical quantities over the entire Brillouin zone. Clus- 
ter extensions of DMFT provide some momentum reso- 
lution 

[Hill, E3j EH E3| and have indicated interes ting 
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momentum-dependent pseudogap behavior [2_ 
I23I |24 , 25l 2fj| . 27 1 . However, most of these works were 
based on studies of two and four-site clusters where the 
momentum resolution is so coarse that information about 
variation around the Fermi surface must be inferred from 
averages involving momentum regions far from the Fermi 




Fermi liquid metal 




FIG. 1: Left panel: Brillouin zone partitioning associated 
with the 8-site cluster and definition of the four inequivalent 
momentum sectors A, B, C and D. The noninteracting Fermi 
surface for t' — — 0.15t and density n = 1 is indicated by the 
gray line. Right panel: sketch of the paramagnetic state DCA 
phase diagram of the Hubbard model on this cluster showing 
a Fermi liquid metal phase, an intermediate phase (hashed 
region) and the insulating state (heavy horizontal line). 



surface. Important exceptions are the work of Macridin 
et al. and Vidhyadhiraja et ai. who studied 16 site 
clusters which have the momentum resolution needed to 
capture the variation around the Fermi surface directly, 
identifying a pseudogap phase at doping 0.05 [28| and a 
critical point defined from the self energy j29[ . We com- 
ment on the relation of our results to these works in the 
conclusion. 

In this paper we use a recently developed continuous- 
time auxiliary field quantum Monte-Carlo algorithm [28] 
to perform direct computations of the momentum vari- 
ation of the electron spectral function and self energy 
around the Fermi surface for the 8-site cluster shown in 
Fig. [1] and for a range of interaction strengths and car- 
rier concentrations. An attractive feature of this cluster 
is that the Fermi surface of the non-interacting problem 
(shown as the heavy line) passes through two symmetry- 
independent sectors (labeled B and C) , so that a coarse- 
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grained approximation to the variation around the Fermi 
surface can be studied directly. We find a two-stage metal 
insulator transition in both the interaction-driven and 
doping driven cases: the Fermi liquid metal and fully 
gapped insulator are separated by an intermediate phase, 
shown as the hashed region in the right panel of Fig. [TJ in 
which sector C is gapped while sector B remains gapless. 
The metal-insulator transition is a momentum-selective 
transition of the kind discussed by Ferrero et al. [27j . 
The one-band Hubbard model in two dimensions reads: 



H 
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with e p — —2t(cos(p x ) + cos(p y )) — 4*' cos(p x ) cos(p y ) 
denoting the electron dispersion and U the on-site re- 
pulsion. We are interested in the low-temperature be- 
havior as a function of interaction strength U and dop- 
ing 6. We use the DCA formulation of cluster DMFT 
m which the Brillouin zone is divided into N = 8 
"patches" defined by the basis functions <fi a (p), which 
are 1 for p in patch a and zero otherwise. The momen- 
tum dependence of the self-energy becomes ~E(p,lu) — » 
^^^(pJEofw), and the frequency dependent func- 
tions £ q (oj), a — 1, . . . , 8 are obtained from the solu- 
tion of an appropriately defined 8-site quantum impu- 
rity model [19(. We solve the impurity model with the 
continuous-time auxiliary field technique (30j with de- 
layed updates [31|. Results are presented for t' = 
(where the less severe sign problem allows us to access 
lower temperatures) and *' = — 0.15*, which is computa- 
tionally accessible but more "generic" than *' = 0, be- 
cause the electron-hole symmetry is broken and the van 
Hove singularity is removed from the half-filled Fermi 
surface. The lowest temperature accessible with the 
computational resources available to us is (3t w 40 at 
U < 8* corresponding (with the conventional high-T c 
band parametrization [32|) to ~ 100K. 

We probe the transition via the sector-specific impurity 
model Green's functions evaluated at the mid-point (3/2 
of the imaginary time interval. This quantity is directly 
measured in our simulations and is related to the value 
of the spectral function of the a sector A a (uj) near the 
Fermi level (co = 0) by (3G a {(3/2) = - J J^A a ( 2 f). 
Figure [2] shows our results for the sectors (B and C) 
which contain the Fermi surface at half filling. In sector 
B (open symbols) we see for U < 6.3* a weakly tempera- 
ture dependent value which is consistent with the nonin- 
teracting Fermi-surface density of states at density n = 1. 
For U > 6.5* the density of states is very small, indicat- 
ing a gapped state. The non-interacting density of states 
in sector C (full symbols) is larger and more strongly 
temperature dependent due to the van Hove singularity, 
which at t' = is at the half filled Fermi level, and at 
*' = —0.15* is about Q.2t below. We see that the temper- 
ature dependence changes sign at U « 5.5* (*' — 0) and 
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FIG. 2: Sector B (open symbols) and C (filled symbols) 
Green's functions, scaled by the number of equivalent sectors 
n a and evaluated at the midpoint of the imaginary time in- 
terval, as a function of U at half filling. The left panel shows 
results for t' — 0; the right panel for t' — — 0.15t. 
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FIG. 3: Imaginary part of the self energy in momentum sec- 
tors B and C for half-filling and t' = 0. The interaction 
strengths are (from top to bottom) U/t = 5, 5.5, 5.75, 6, 
6.25, 6.375, 6.5, and 6.625. 



U k 5.8t (t! = — 0.15i); for larger U the evolution is to- 
wards an insulating state. Thus sector C undergoes a gap 
opening transition at a smaller U than sector B, both for 
t' = and t' = — 0.15i. We emphasize that this two-stage 
transition is directly visible in our data and is a result en- 
abled by the momentum resolution of the 8-site cluster. 
It is different from the result of single-site DMFT, which 
predicts a single transition at a large U ~ 1.5W with W 
(here 8t) the full bandwidth and from the result of 4-site 
cluster calculations, which predict a single transition at 
U - 4.2* (DCA 26]) or - 5.5* (CDMFT (H). 

Figure [3] displays the imaginary part of the Matsubara 
axis self energy in sectors B and C for several U in a 
range containing the two transitions. We present results 
obtained for t' = both because the absence of the sign 
problem allows us to obtain more accurate data and be- 
cause the particle-hole symmetry implies ReE = which 
simplifies the discussion, but our results for *' = —0.15* 
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FIG. 4: Main panel: B and C sector dopings (summed over 
spin and sector degeneracy) and total doping as a function of 
chemical potential for U/t — 7, t' — — 0.15t at fit = 20. Inset: 
expansion of electron-doped metal insulator transition region 
showing the temperature dependence. 



are consistent. The smaller U curves of Fig. [3] show the 
low frequency behavior expected in a Fermi liquid phase 
Im£(iu; n ) ~ — w„(Z _1 — 1) with quasi-particle residue 
Z decreasing as U is increased towards the transition. 
The larger U curves reveal the low frequency behavior 
ImE(iw n ) ~ —A 2 /uj n expected in an insulating phase, 
with A decreasing as U is decreased towards the metal 
insulator transition value. In the intermediate regime 
5.5 < U/t < 6, Sc exhibits an insulating behavior while 
Us appears consistent with the Fermi liquid expectation. 

Figure [?] presents the doping driven transitions for 
U = 7t, large enough that at half filling all sectors are 
gapped. We denote by \l the chemical potential minus 
the Hartree shift Un/2. For t 1 = — 0.15t hole-doping 
leads to a wide range of chemical potentials where sector 
C remains insulating while sector B is doped. On the 
electron-doped side, the two transitions occur almost si- 
multaneously and the onset of doping in sector B appears 
very sudden. However, the inset to Fig. 2] shows that 
even on the electron-doped side there is a small window 
of chemical potential (// ~ 0.7t) where sector C develops 
a gap but sector B remains gapless. In the model with 
t' = (not shown) the £?-sector transition (at (3t = 40) 
occurs near \fi\/t = 0.5 and the C-sector transition near 

|mIA = i. 

The particle-hole asymmetry in the model with t' = 
— 0. 15t is evident also in the self energies. We find that on 
the electron-doped side, at high dopings S > 0.15 the two 
sectors have the same quasi-particle residue (Z a « 0.5 
for S = 0.23); as doping is decreased, Z decreases in both 
sectors, but at a slightly faster rate in sector C. Down 
to a doping of 5 m 0.08, Z is reasonably well defined 
in both sectors, is larger than 0.2 and is only weakly 
anisotropic. For lower dopings in the vicinity of the C- 
sector gapping the behavior is not Fermi liquid like and 
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FIG. 5: Imaginary part of B-sector self energies for the doped 
Hubbard model at U = 7t, t' = -0.15t and (it = 20. Traces 
indicated by open circles, triangles and diamonds correspond 
to jj, where both sector B and sector C are gapped; traces 
marked by full circles, triangles and diamonds correspond to fj, 
where sector C remains undoped and sector B is doped, while 
traces marked by asterisks, crosses and plus signs correspond 
to fi where both sectors are doped. 



the C-sector self energy becomes much larger. Figure [5] 
shows that the hole-doped side exhibits a richer behav- 
ior. An analysis of the real and imaginary parts of 1/E 
(not shown) reveals that at — OAt < fi where the model 
is insulating (gaps in all sectors) the low frequency self 
energy is dominated by a pole E Q ~ A^/(iuj n — uu^) 
with pole strengths Ag/t w 3 and Ac/t « 5 and po- 
sition ujp shifting with chemical potential. This shift 
is responsible for the changes evident in the three un- 
doped traces. In the sector selective regime where B is 
doped and C is not, the B self energy is that of a 'bad 
metal' (lim^^Q lmS(iu) n ) =/= 0): any Fermi liquid coher- 
ence scale would be far below the temperature regime 
we can access. However, as the doping is increased into 
the regime where all sectors are ungapped the self energy 
evolves towards Fermi liquid behavior (with Z a w 0.4 for 
S = -0.23). 

The results are intriguingly similar to data on high- 
T c materials where a growing body of evidence indicates 
a metal-insulator transition in the hole doped materi- 
als apparently occurring by a reduction in the size of 
the fc-space regions which can support mobile particles, 
while the electron-doped materials behave more conven- 
tionally (apart from the effects of long ranged commen- 
surate density wave order not included here). We be- 
lieve that our coarse momentum resolution approximates 
the continuous reduction of Fermi surface area by a two 
step transition. Preliminary results for a 16-site cluster 
with Fermi surface crossing 4 inequivalent sectors indi- 
cate a 4 stage metal-insulator transition consistent with 
the gradual opening of a gap near (0,7r). The electron- 
hole asymmetry suggests an important role played by 
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the van Hove singularity in the sector-selective transi- 
tion. On the electron-doped side the sector-selective re- 
gion shrinks with increasing particle-hole asymmetry and 
may completely vanish at t' = — 0.3t. 

One can view the differentiation between sectors B and 
C as an orbitally selective transition of the effective im- 
purity model 27J. ft is known that orbitally selective 
transitions may be destabilized by a Kondo coupling to 
other sectors; however in our case the gapped sector C is 
two-fold degenerate and we believe that the two electrons 
in this sector are strongly antiferromagnetically coupled 
(or form singlets) so that the state is stable against a 
Kondo coupling. We therefore expect that the effect sur- 
vives beyond the DCA approximation we have used. 
Our results are in general agreement with previous 
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cluster DMFT work 20, 28] indicating a 'pseudogap' phe- 
nomenon more pronounced for hole-doped than electron- 
doped materials. Our improved momentum resolution 
confirms that the effects inferred from two and four-site 
calculations 0, ED, Q HJ [H, ^ indeed reflect be- 
havior associated with near-Fermi-surface states, while 
the improved temperature resolution and the examina- 
tion of a range of interactions and dopings reveals the 
sector-specific metal-insulator transitions which underly 
the pseudogap observed by Macridin et al. [28]. Note 
that while Ref. [lU reported a large density of states 
near (0, it) for 5% doping both for t' = and t' = — 0.3t, 
we find (for t' = and t' = — 0.15t) that dopants first ap- 
pear near the nodes (sector B) with sector C remaining 
gapped. However, after the initial doping, the sector C 
occupancy increases rapidly, so the difference may be a 
doping- level effect. The self energy transition studied in 
Ref. [29( can be seen in Fig. [5] and the corresponding data 
for sector C (not shown) as a change from IvaT,{iuj n ) ~ uj n 
to a flatter w-dependence as 8 drops below ~ 0.15. This 
change occurs well before the sector C gapping transi- 
tion (5 ~ 0.08) and its connection to the sector-selective 
transitions we find remains to be clarified. 

The important open questions concern the physical ori- 
gin of the gapping behavior. We find cluster antifcrro- 
magnetic correlations which are large and grow rapidly in 
the vicinity of the metal-insulator transitions. However, 
in the model with t' = — Q.15t the correlation functions do 
not change dramatically across the sector-selective tran- 
sitions, suggesting that the phenomenon may not be an- 
tiferromagnetically driven, perhaps supporting the inter- 
pretation of Lauchli et al. [12j . Optimized implementa- 
tions of the hybridization expansion technique [33| should 
permit examination of the cluster cigcnstates along the 
lines of Ref. [261 ] , and thus provide direct information on 
the role of singlet formation and antiferromagnetism. 
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